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Abstract 

We investigate the fluctuation contribution to the specific heat of a two- 
dimensional superconductor with a non-Fermi normal state described by a 
Anderson Green's function C7(k, iuj) = u~ a /{iu> — £k) 1_ °- The specific heat 
corrections contain a term proportional to (T 2a — T 2a )~ 1 and another logarith- 
mic one. We defined a coherence length as function of the non-Fermi parameter 
a, which showed that a crossover study between BCS and Bose-Einstein con- 
densation is possible by varying a in an interval -j- a cr . By comparing our 
theoretical results with the experimental data for HTSC materials, we reob- 
tained the value for a, corresponding to such systems, of the order 0.3 -f- 0.45. 
We also reobtained the critical temperature for such a superconductor using 
the Thouless criterion. 
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INTRODUCTION 



After the discovery of high critical temperature superconductors (HTSC) a great number 
of experimental works showed that the usual Fermi liquid description is unappropriated for 
these materials. As a consequence there are several phenomenological mo delM proposed in 
order to explain the high critical temperature T c and the nonmetallic behavior of such mate- 
rials. An important role in the mechanism of superconductivity is related to the fluctuation 
effects which are very important in the critical region. One of the parameter that characterize 
the normal state properties of a HTSC is the specific heat. Thompson and Kresini, showed 
that in the absence of a magnetic field B = 0, in the critical region (T — > T C ,T > T c ) the 
corrections in the specific heat obtained due to the fluctuation effects are proportional to 
(T — X^ 1 / 2 . These corrections correspond to the classical region of fluctuation!'!. Recently 
Grosui reconsidered the problem for a two-dimensional superconductor with a logarithmic 
density of states and he showed that the corrections due to the fluctuation effects in the 
specific heat consist in two terms, one proportional to (T — T c ) -1 and another logarithmic 



one. 



In a previous papeiH we derived the form for the specific heat jump in a non- Fermi 
superconductor described by the Anderson model. In our calculations we considered that 
the normal state Green's function has a similar form as the one from ID interacting fermion 
systems^. A similar choice for the Green's function G(k, iu) = uo~ a /{iuo — ^k) 1 ^ where u c 
is a frequency cutoff introduced to makes the Green's function dimensionally correct and 
a is a parameter related to the anomalous Fermi surface, was made by Chakravarty and 
Anderson0 in order to study the interlayer mechanism of the cuprate-oxides. The model has 
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been developed by SudboS who showed that the parameter a has to satisfy < a < 0.5. 

The aim of this work is to calculate the specific heat corrections of a two-dimensional 
non-Fermi system using the Anderson model. In Section II we derive the particle-particle 
polarization for the model and we reobtain the critical temperature T c . We introduce using 
the fluctuation propagator a Ginzburg-Landau coherence length extrapolated to zero tem- 
perature, which in our calculations is function of the parameter a. In Section III we evaluate 
the correction terms in the specific heat due to the fluctuation effects. In Section IV we 
analyze our results. 

PARTICLE-PARTICLE POLARIZATION 

In order to evaluate the corrections in the specific heat inside the critical region above 
T c first we have to calculate the Green's function connected with the Cooper pairs. It was 
showed by SchmidEl that this Green's function is given by 



^ + n( q ,^ 



where V is the interaction strength and II(q, u n ) is the usual particle-particle polarization 
given by 

r d 2 k 

II(q, u m ) = TJ2 J G ( k 7 ^n)G(q - k, uj m - tu n ) (2) 

with 

G(k,u n )= (3) 

The behavior of the Green function given by Eq. (Q) for uo — > has been analyzed in Ref. [12 
and the anomalous a;-dependence of it in the interval — u c , uo c gives a softer singularity than 



1/u specific for the Green function of the Fermi liquid. 

We remark that in the usual Fermi liquid theory (a = 0), for a constant density of states, 
in a three dimensional system we reobtained the standard result (See Ref. [13]) 



D(q,u n ) = — — ■ — j (4) 

N [t + a\uj n \ +£VJ 

where N Q is the density of states at the Fermi level, r = (T — T c )/T c , a = ir/8T c and £ is the 
Ginzburg-Landau coherence length extrapolated to the zero temperature. 

In order to calculate the particle-particle polarization we will rewrite Eq. (Q) in the 
following form 

f d 2 k 

II(q,^ m ) = J ^-p5(k,q,w m ) (5) 

where 

LJ~ 2a 

The evaluation of the sum over the Matsubara frequencies in Eq. (|J) is made using a contour 
integral which involves the Fermi-Dirac distribution function, with the specification that in 
this integral we have to deal with two branch cuts related to the complex function involved 
in the summation. Another observation which has to be made is the fact that a correct 
evaluation of this sum is possible only on the assumption that we work at low frequencies 
ui m and small wavevectors q. In this way we get 

u~ 2a sin [ir(l — a)} 



S(k,q, io n 



TV 



x {-(2e ]l ) 2a - 1 B(l - 2a, a) 

+ fI ( -' r ^ ^ M, ' + " 1 



+ (a - l)(v F qcos9 - u m ) -(2e k ) ia " 2 J B(2 - 2a, a - 1 



+ ^£<-u 



,«-l/2 



m=0 



+ 1^-3/2^ a ~ 3 / 2 



K a - 3 /2[£kP(m+ 1)] 



(a — l) 2 (t>_Fgcos#) ; 



+ 2r(a-i) 



7T 



m=0 



-(2£ k ) 2Q ^B(3-2a,a-2) 

(m+l^-V2 ^-3/^M( m + 1 )] 
(2e ) Q ~ 3 / 2 



(7) 



where for the validity of this formula a should be between and 0.5, B(x,y) is the beta 
function, T(x) is the Euler's gamma function and K u (z) is the Bessel function of imaginary 
argument. 

Using Eqs. (||) and (0) is very simple to evaluate the particle-particle polarization 



Il(q,cu) 



2jV sin [tt(1 -a)} 



7T 



2— B(l - 2a, a) + ^ 1 - 2^ /T 

W 



+ 



— a) 



2 2a - 2 5(2 - 2a, a - 1) ^ D ^ 

UJ r 



2a 



2a - 1 



r(a-l)rfr-l/2)l-2^ ^ 



2 2-2 C 



2a-l' 



+ 



(f F g) 2 (l - a) 
4w 2 



-C(2a-1) - 

2 2a - 3 5(3-2a,a-2) /u; D \ 
2a - 2 



2a-2 



+ 



'2r(a-2)r(a- 1/2) 



T 2 fa- 1^ 



7T 



1 - 2 3 ~ 2a 
2 3 " 2a 



2a-2 
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where ((x) is the Riemann's zeta function. 

In the case of zero frequency and zero wavevector n(q = 0, iu — 0) we can obtain the 
critical temperature for normal state-superconductivity phase transition using the Thouless 
criterion 



1 + VReU{q = 0, iu = 0) = 



(9) 



It is easy to see that we reobtained the previous result from Ref. |T4] . |T5 



T 



2<> 



C{a) 



D(a)u%- C 



2« 



A(a)X 



(10) 



where A = N V, C{a) = T 2 {a){l - 2 l ~ 2a )({2a), A{a) = 2 2a sin [tt(1 - a)]/ir and D(a) = 
B(l - 2a, a) /2a. 

Using Eqs. (§) and ( fT0|) the fluctuation propagator can be written in the following form 



D-\q,Lj) = N A{a)C{a 
+ N A(a 



rp ^ 2a i rp x 2d 



LO, 



UJ, 



-M(a)(l-a)+ 1 



T 



4T 2 



iV(a)(l-a) 2 (^) 



where 



M(a) 



+ 



_ 5(2 -2a, a- 1) /c^X 2 "- 1 T, 
2(2a - 1) W / _ 
r(a-l)r(a- 1/2) 



2v^F 



1 _ 2 z - 2a )((2a - 1) ( — 



2a 



'2r(a - 2)r(a - 1/2) 



2a-2 



+ 



4(2a - 2) V cj c / 
„ _ 1 /o^ 

+ r 2 («-r 



7T 



L ^—C(2 - 2a) (i£ 
4 W 



2a-2 



which can be written also in the following form 



(12) 



(13) 



£>(q,u;) 



N [br + iauj + £ 2 {a)q 2 l 



with a = (1 - a)M(a)A(a)/T, b S A(a)C(a)(T c /a; c ) 2 , r = (T/T c ) 2a - 1 and 



,2/ x ~ ^(!-«) 2 

^ = 4„ 2 



iV(a)A(a) 



(14) 



(15) 



If we compare Eqs. fli~4 ) and ([|) we can define a Ginzburg-Landau coherence length 
extrapolated to zero temperature, which now is a function of the non-Fermi parameter a. 
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As it is known an important parameter related to the crossover theory between a BCS like 



of such a parameter is about 10 3 -^10 4 and in the Bose-Enstein condensation limit is zero. 
Experimental data showed that in the case of HTSC materials this parameter is between this 
two limits, for example in the case of YBa2Cu307_<5 is about 5-hlO. In Fig. ([[]) we plot the 
variation of the parameter kp^ Q (with the critical temperature corresponding to the optimal 
doping of YBa2Cu307_,5) as a function of the non-Fermi parameter a. The limit a — > 
correspond to the overdoped limit where the usual Fermi liquid theory still works for the 
normal state properties. The other limit a — *> a cr correspond to the underdoped limit where 
the usual Fermi liquid theory breaks down. In this limit as we can see from Fig. ([!]) the 
coherence length is going to zero, which shows that here a Bose-Einstein condensation of 
strong coupling pairs (£ — > 0) is more close to the real picture of the system. The intermediate 
limit (a — > 10), which correspond to the optimal doping of YBa2Cu307_,5, shows that the 
non-Fermi parameter a should be around 0.30-^0.45, result which is in agreement with our 
previous calculation of the specific heat jump at the critical temperature!. 



The specific heat corrections due to the fluctuations around the critical temperature, 
will be calculated by having an additional term in the thermodynamic potential. Inside the 
critical region only the zero frequency term is important 



formalism and a Bose-Einstein condensation formalisirJ-3 is /c_f£o- in the BCS limit the value 




FLUCTUATION CONTRIBUTION TO THE SPECIFIC HEAT 



Q /z =T£ln [D- 1 (q,io; = 0) 



(16) 



q 
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Using the well known thermodynamic relations S = —dQfi/dT and C = TdS/dT the specific 
heat correction close to T c is 

^d^fifdrV dn fl d 2 r 
Cfl -~ Tc ^{dT) ~ Tc ^TdT 2 (17) 

From Eqs. flU) and (0) we get 

ba 2 T 2a ba(l-2a). ( 2mE F £ 2 (a) T 2a \ 



Cfl = ^T-^TT + „ ZJTi ^ 



TiQ(a)T 2a -T 2a 2vre 2 (a) \ b T 2a - T 2a J v ' 

The result contained in Eq. ( |i8| ) is valid only in the critical region when we consider the 
classical fluctuations (region where the random phase approximation is still valid). Close to 
the transition point the interaction between electrons and Cooper pairs are very important, 
so we expect that in some way the electronic Green's function is affected. As a consequence 
there is a contribution 511 to the particle-particle polarization II in Eq. @, which can be 
written as 

5IL = 2T 2 W 7^G 2 (k, tu n )G 2 (-k, -co n ) f 7^^(q, iu = 0) (19) 

If we perform the integrals over k and q, and if we sum over the Matsubara frequencies u n 
we get 

With this correction the particle-particle polarization becomes lb = II + 511 and the fluctu- 
ation propagator can be written as 

D; 1 = IT 1 + 511 (21) 

Using this definition very close to the transition temperature the significant contribution due 
to the fluctuations will be 
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n t ^ g In 



cln 



2mE F e {a) 



br 



(22) 



where g is a constant which contains the critical temperature T c and a cutoff wavevector of 
order of q F (the Fermi wavevector) and c is the constant in front of the logarithmic function 
in Eq. (ED). With this thermodynamic potential the fluctuation contribution to the specific 
heat close to the critical temperature becomes 



a 



4a 2 g 



br 



i -i 



2a(2a - l)g 



t In 



2mE F Q {a) 
br 



i -i 



(23) 



As we can see from Eqs. ( p23[ ) and (pif) the effect of the interaction between electrons and 
fluctuations should be very important, at least close enough to T c . The only question is if 
the formula from Eq. (|23"D is valid in the framework of the performed approximations. There 
are several observations which make questionable this formula. The most important one is 
the fact that in Eq. flT9| ) we used instead of a renormalized fluctuation propagator the bare 
one, which gives us a different behavior as a function of the wavevector q for SU. In fact 
Eq. ( |T9"D should be a self consistent equation for SU, which in our work is solved in the zero 
order. 



DISCUSSION 

We evaluated the contribution of the fluctuation contributions in the specific heat of a 
non-Fermi system described by the Anderson mo del We derived the form of the particle- 
particle polarization, and we showed that the critical temperature can be reobtained using 
the Thouless criterion. Using the fluctuation propagator D(q, uo) we introduced a coherence 
length £,o(a), which inn our calculations is function of the non- Fermi parameter a. By varying 
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a between zero and a cr we showed in the case of YBa2Cu30y_5 that the system can be 
described using the BCS theory or by the Bose-Einstein condensation of preformed pairs. In 
the case of optimal doping we reobtained for a the same values as in our previous paper! 



a = 0.3^0.45. As a consequence we consider that the numerical analysis in Ref. P] , |15l should 
be made on a scale a = 0.3 -v- 0.45, not on the scale very close to the value a = 0. 

In the third section of our paper we calculated the fluctuation contribution in the specific 
heat close to the critical temperature T c , related to the classical fluctuation region, where 
the random phase approximation still remains valid. The result obtained by considering the 
fluctuation contributions give two different terms, the significant one being proportional to 
(T 2a — T c 2a ) _1 . A similar result can be obtained if we consider the critical field H C 2(a,T) 
obtained by Muthukumar et al.0 using the Ginzburg-Landau theory. Close to the critical 
temperature T c , where interaction between electrons and Cooper pairs becomes very im- 
portant we calculated the behavior of the specific heat by introducing a correction in the 
particle-particle polarization. The corrections obtained due to this additional term is very 
important, the divergence in Eq. ( PB"D being more powerfull than the one in Eq. flT8|). All 
these results are valid in the approximation of the "box corrections" for the fluctuations 
propagator. 

Finally we mention that the superconducting state developed in Ref. [1^ using this non- 



Fermi model presents some differences to the results from Ref. 11 because of the phase factor 



in the Green function. However, using the model proposed in Ref. 11 in the limit a — > it 



was 



showedS that the critical temperature is identical to the BCS expression. 
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